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abstract 


In this thesis, a 6T0 converter having six GTOs is 
considered for state space model development . This model has 
been developed, based on the graph theoretic approach which is 
conceptually simpler and leads to efficient formulation of the 
converter equations corresponding to al 1 possible modes of 
operation. The state equations are derived using network cutset 
matrix which can be modified conveniently depending upon the 
converter conduction pattern. The results of steady state test 
simulation are presented for both modes of the GTO converter 
operation (rectifier / inverter ) to i 1 lustrate the validity of 
model and capability of program developed. 

The performance of the HVDC link with GTO inverter has 
been studied under both normal and abnormal condi t ions .The system 
fault studies are change in current reference and ac voltage 
dips at both end of the converter . The harmonic analysis 
are performed for both normal and abnormal conditions. 



CHAPTER 1 

INTRODUCTION 

High voltage direct current (HVDC) transmission 
has become a viable alternative to alternating current 
transmission for bulk power exchange. The earlier HVDC 
installations used mercury arc valves for power conversion. 
After the advent of thyristor, the most of the HVDC installations 
employ thyristor valves in converters which operate under 
natural (line) commutation. The problems of reactive power 
requirement and harmonic generation in case of line commuted 
converters are well documented in the literature C13. These 
problems assume greater significance in case of the dc 
transmission and particularly when the power transmission 

capability of the line is being continuously upgraded. 

Traditionally, the harmonics generated by converters 

( rectif ier/invereter ) are not allowed to enter the associated ac 
system by using tuned filters at the converter ac bus. The J 
filters can effectively shunt only the characteristic (6n+1 ! 

,n=1,2,«...) harmonics while the uncharacteristic harmonics may 

1 

still enter the ac system and may cause harmonic instability £211. 
The reactive power requirement of converters which may be j 
typically 50% of the real power under normal operation, is met to a : 
certain extent by installing shunt capacitors. Also, the harmonic ' 
filter which produces reactive power at fundamental frequency, ; 
provides some reactive power support. To circumvent the 



problem of reactive power requirement and harmonic generation, 
the use of farced commutation in thyristor converters has 
been envisaged particularly with reference to inverter 
operation C33. Pulse Width Modulated (PWM) ac-dc thyristor 
converters have been suggested to overcome the drawbacks 
of reactive power and hormonic generation. Although, it 
improves power factor and reduces harmonics, but the 
requirement of large number of auxiliary power devices and LC 
elements for the purpose of forced commutation of thyristors make 
these converters heavy, expensive and less reliable. Recently 
the availability of Gate Turn Off (GTO) thyristors with improved 
reliability and high ratings have made them attractive to be used 
in PWM ac— dc converters. With this, better power factor at ac 
side, less distortion of the line current and reduced ripple 
in dc current are achieved. As a GTO requires no commutation 
arrangement in the power circuit and can be turned off by a 
negative gate signal, its use as a power device, in place 
of a thyristor with forced commutation arrangement, greatly 
simplifies the converter configuration .The application of gate 
turn off(GTO) thyristor in low power circuits has been reported in 
literature C4 3. It may be interesting to examine its application 
in HVDC converters ( parti cularly inverters > using pulse width 
modulation (PWM) scheme. 

1.1 GTO CONVERTER WITH PWM 

In view of the advantages, a number of PWM schemes for 
ac-dc converters have been listed C4,5 3. All these schemes have 
their own merits and demerits as regards to performance on source 



and load side. The comparative study of various PWM schemes have 
been made for three phase PWM ac-dc converter employing GTO 
thyristors connected to an active load C43. The basis for 
comparison is the source side performance such as power 
factor, distortion of line current, percentage of lower order 
harmonics in line current with fundamental component and 
the load side performance such as ripple content in the output 
current. Mainly five different types of PWM schemes, vizs(i)Equal 
(ii) Triangular (iii) Stepped (iv) Sinusoidal and (v) Inverted 
sinusoidal are in vogue. Among these schemes. Equal Pulse Width 
Modulation (EPWM) scheme has been found to offer the best overall 
performance C43. 

The various pulse width modulation schemes can be used 
with converter configuration which employ either all 6 GTOs or 3 
GTOs and 3 thyristors. While the latter offers economy in cost as 
GTO is a very expensive device compared to a converter grade 
thyristor but it requires large number of switching operations 
per cycle thus increasing the switching losses. Furthermore in 
the latter configuration firing pulses are to be reset for a 
different gating sequence whenever the operation has to be changed 
from rectification to inversion mode or vice-versa. Since turn on 
and turn off of a GTO is under an external control, the gating 

i 

patterns (which turn on & turn off a GTO) could either be 
symmetrical or unsymmetrical C43. In view of the above 
description, the six GTOs converter circuit employing symmetrical! 

ft 

gating sequence based on equal pulse width modulation scheme, has 
been considered in this thesis for various investigations. 



1.2 OBJECTIVES OF THE THESIS 


The primary objectives of this thesis are s 

(i) to examine the feasibility of using GTO converters 
as inverter in HVDC transmission system. This investigation is 
carried out through detailed dynamic digital simulation of a two 
terminal HVDC link. 

(ii) development of a suitable mathematical model of a GTO 
converter for the purpose of simulation. 

1.3 OUTLINE OF THE THESIS 

In chapter 2 , a brief review of the literature on the 
PWM converters and various simulation techniques is presented. 

The detailed modelling of GTO converter is given in 
chapter 3. This chapter also includes firing schemes 
(i.e. symmetrical and unsymmetri cal ) and effects of source 
inductance. The state space model is developed for the GTO 
converter with active load and RC transient suppressor 

circuit. Using the concept of graph theory, a generalised 

set of differential equations is formulated. The results ol 
simulation are also presented in this chapter for botl 
rectifier and inverter operation. 

Chapter 4 describes the modelling of the constituent^ 
of an HVDC system in brief. Emphasis is laid on interfacing of GT( 
inverter and its controls in place of thyristor inverter usin; 
appropriate interfacing variables. The model developed in chapte : 
3 for GTD converter is used to study the performance of HVDC link; 
A computer program is developed to demonstrate the feasibility o 
GTD inverter in HVDC link. The disturbances considered, to stud 



the dynamic behaviour, the change in current reference setting and 
ac voltage dips at both rectifier and inverter terminals. 

The conclusions and future scope of the work are 
in the chapter 5. 


given 



CHAPTER 2 

LITERATURE REVIEW ON GTO CONVERTERS AND RELATED TOPICS 


A brief review on pulse width modulated GTO converters 
and simulation techniques are given in this chapter. 

£.1 PULSE WIDTH MODULATION 

Line commuted converters are widely used in 

power industries because they are simple and reliable. 
However, these converters deteriorate the power factor at the 
source specially at large phase angles and they also 
introduce lower order harmonics in the source current which can 
be filtered out only by a large size LC filters. The HVDC 

transmission, using line commuted thyristor converter 

suffers from the problems of lagging reactive power 
requirement at the inverter and the commutation failure 

under abnormal conditions. Pulse Width Modulated (PWM) 

ac— dc thyristor converters have been suggested to overcome 

these drawbacks. Although, they improve power factor and 

reduce hormonics but require large number of auxiliary power 
devices which make the converter heavy, expensive and less 
reliable. Availability of high power GTO thyristor which can be 
turned on and off at desired instant, offers a suitable 
alternative. The PWM techniques employing forced 

commutation can make the fundamental power factor unity or even 
leading in the entire control range of the output voltage. But 
few voltage pulses per cycle introduce more ripples in the 



output current and lower order harmonics in the source current. 
The width of output voltage pulses is continuously modulated 
throughout the ac cycle according to the modulation of a 
modulating signal and the PWM scheme is named after the 
corresponding modulating signal. A number of PWM schemes 
such as Equal Pulse— Width modulation (EPWM) ,Trangular 
Pulse-Width Modulation (TPWM) ..Stepped Pulse-Width Modulation 
(STPWM) , Sinus iodal Pulse-Width Modulation (SPWM) and Inverted 
Pulse— Width Modulation (ISPWM) are reported in the literature 
143. 

Z.Z CQMPERATIVE PERFORMANCE OF VARIOUS PWM SCHEMES 

The detailed comparative study of source and load 
side performances of a three phase PWM ac— dc converter 
employing different PWM schemes is given by Khan C43. The effect 
of varying the number of output voltage pulses is also given. The 
switching instants are determined by comparing a trangular 
carrier signal Vca with a modulating signal Vm as shown ir 
Figure (2.1). The source is connected to the load for the duratior 
when Vm>Vca. The frequency of carrier signal depends on tht 
number of pulses(P) required in one cycle of ac source. Th* 
output voltage can be varied by varying the peak value, Vmm o( 
the modulating signal with respect to the peak value, Vcam o-f 
carrier signal. The ratio Vmm/Vcam is termed as modulatioi 
index <m). The expression for different modulating signals ani 
the corresponding firing angles (<x*s> and extinction angles 
for Kth pulse can be given as follows C4I 


Vca=+VcamC ( P /n ) * ( wt ) - ( 2K-1 ) 3 




TPWM 

FIG.2.1. different modulation schemes 



Where + stands for positive slope and - stands for negative 
slope. The expression for firing and extinction angles for 
different modulation schemes are given in Table 2.1. It is 
reported in C4U that EPWM is superior to the other PWM 
schemes because of linearity between dc voltage ratio ( y ) and 
modulation index (m) which is advantageous when the converter is 
used in a closed loop schemes. It is also reported that EPWM 
scheme results in least current ripples among various modulating 
schemes both in rectification and inversion operations. EPWM has 
only nP+1 and nP-1 source current harmonics. EPWM with 12 pulses 
per cycle is used in this thesis. 

2.3 GftTE TURN-OFF THYRISTOR (GTO) 

GTOs were offered at the start of sixties with 
higher voltage capability but lower current capability. But since 
last few years they are available with ratings as high as 4.5 
kV,3000 Amp. C63. GTO thyristor is a silicon semiconductor power 
switch similar to an ordinary thyristor, but with a difference 
that it can be turned off (as well as on) by appropriate pulse at 
the gate terminal. Unlike an ordinary thyristor, which needs 
a higher power bulky circuit for turn off, GTOs can be turned 
off from a low power electronic circuit. As may be 
expected , GTO thyristor brings with its some drawbacks as i 
well as advantages. Since turn off time of GTO thyristors isj 
low, therefore the rate of change of anode current is very high. 

I 

Unless the stray inductance is kept to a minimum, an excessive 
voltage will be developed across the inductance. The layout of 
snubber circuit requires special circuit to 


minimise stray 



s Expressions for different modulating signals and corresponding firing/extinction 
angles for Kth pulse. 

Modulating signal v m Firing angle <x R T Extinction angle f v 



iductances £7,8,91. There are two basic types of GTO 
hyristars, the anode shorted emitter and the gold diffused 
evice. Due to some advantages, the gold diffused devices are 
refered USD. Both GTO thyristors and conventional 

hyristors are latching devices. Once the anode current reaches 
o a certain latch— on value , the device will automatically go 
nto a full conduction state. Only difference is that 
onvent ional thyristors can be brought back to the blocking 
tate only by reducing the forward anode current below the 
lOlding current and GTO thyristor can be brought back to the 
blocking state by applying a negative gate signal £6,101. 

>.4 PWM ftC-DC GTO CONVERTERS 

Basically, there are two types of PWM ac— dc converters 
is given below. 

1.4.1 Three GTOs and Three Thyristors PMM AC-DC Converter 

Three GTOs and three SCR thyristors converter 
configuration is proposed by INABA et al £111 which is given in 
Figure (2.2). Since GTOs are more costly than the ordinary 
thyristors, this converter has lower cost than the converter 
using six GTOs. The main problem is that a modification of 
gating sequence will be required to operate the converter in 
inversion mode. Overlapping between the conduction period of 
the outgoing and incoming GTOs is also provided to ensure the 
continuity of load current. 

5.4.2 Six GTOs PMM AC-DC Converter 

The advantage of using six GTOs is that one can use 
optimal PWM control schemes, where particular harmonics can be 
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FIG. 2. 2. CONVERTER CIRCUIT GIVEN BY 
INABA et al. 




inimized, but it increases the switching times leading to 
ncrease in switching losses. There is no requirement of 
odification of gating sequence to shift operation from 
edification mode to inversion mode or vice-versa. Because 
he GTO is expensive compared to thyristor , the six GTOs PWM 
ic— dc converters are more costly than three GTOs and three 
thyristors convertors. The six GTOs configuration is given by 
(atakoa et al C 12U which is given in Figure (2.3). 

Z.5 SUPPRESSION OF SWITCHING SPIKES 

When the source current is forced to zero, the 
source inductance produces voltage spikes which must be kept 
to the acceptable level. An RC filter is often connected on ac 
side for this purpose. This leads to substantial wastage 
of power, consequently PWM converters have lower efficiency. 
Viriya et al C133 have suggested an alternative circuit for 
suppresion of the voltage spikes, which improves the 

converter efficiency significantly and is -given in Figure 

(2.4). An auxiliary circuit capacitor is used to clamp the 
commutating voltage spikes at the desired level. The reactive 
energy stored in the source reactance, at the time of commutation 
of each GTO in the main circuit is passed through to the 
capacitor in the auxiliary circuit. The energy stored in the 
capacitor is then discarged in pulses to the dc load during on 
state of the GTOs in the auxiliary circuit. By this method the 
snubber capacitor becomes smaller or even unnecessary. 

2.6 SIMULATION TECHNIQUES 

For study purpose, there are two type of models 





namely, physical and analytical. Physical model has 
advantage of experimentation but it is less flexible. The 
mathematical description of the system characteristic in 

the form of algebraic and differential equations has the 
advantages such as flexibility for analysing with different 
input— output parameters, usefulness for studying various system 
alternatives etc. Although some assumptions have to be made for 
analytical model. There are four basic simulation techniques 

which have been employed for the study of power semicoductor 
circuit C53. 

(1) Digital Simulation 

(2) Analog Simulation 

(3) Breadboard Simulation 

(4) Parity Simulation 
2.7 DIGITAL SIMULATION 


For digital 

simulation, the 

system under 

study 

is characterised by a mathematical 

model . 

The dynamics 

of the 

system can be discribed 

in the 

general 

state space 

form as 

X=CAIX+CBIU 


(2.1) 



Y=CCDX+I:D3U 


(2.2) 




where X is the set of the system state variables, U is the 
vector of forcing functions, Y is the set of output variables of 
interest. The matrices CAI,CBI,IICI] and CD3 are system matrices and 
depend upon physical characteristics and parameters of the 
system. The state equation (2.1) is solved using numerical 
integration and set of output variables are calculated at each 
time step by equation (2.2). The varying circuit configuration 



leads to the time 

dependent nature of 

the system 

matrices 

and 

thus formulation of 

the state 

and 

output equations are 

not 

straight forward. 

The dynamic 

representation 

of the power 

circuits has been 

reported in 

the 

literature 

C143 based 

on 


(a) constant topology approach 

(b) varying topology approach 


2.7. 1. Constant Topology Approach 



In 

constant topology 

approach the 

network 

structure 

i s made 

time 

invariant 

for 

all modes of 

operation by 

suitable 

representations 

of the devices. 

The 

switching 

devices 

are 

represented 

by a 

low resistance 

when 

conducting 

and by 

high 

resistance 

when 

not conducting. 

This 

model runs 


the risk of numerical instability because of widely varying 
time constants if the step size of integration is not suitably 
selected. The problem relating to varying time constant is 
solved with a new model for switching devices. In the 
nonconducting state of the switching device the reverse 
resistance is arranged in series with an appropriate value of 
fictitious inductance so that this time constant is larger 
than the real time constant in the circuit. The switching 
devices may be represented by a combination of forward and 
backward resistance to that structure of the signal flow graph is 
independent of the mode of operation. Another way to represent 
the switching devices is to adjust voltage source. 

2. 7. 2. Varying Topology Approach 

In this approach the switching devices are so modeller 
that the network exhibits a time varying structure. Th« 



equations are formulated from the reduced network structure for 
the particular mode of operation with time varying structure. 
Spurious time constants are not introduced in mathematical model 
and also the order of the system matrices is minimum for a 
particular mode of operation. Various models reported in 
the literature represent the switching device during off 
state by an infinite impedance (i.e. zero current through 
device) and during on state either by zero impedance ti.e. zero 
voltage drop across deivce) or by small resistance. 
Sachchidanand et al C153 have represented the converter based on 
the graph theoratic framework for the detailed dynamic 
simulation of a three phase bridge converter. This approach not 
only formulates the system state equations efficiently for all 
possible modes but also gives the expression of dependent 
variables. This enables to analyse the converter during both 
normal and abnormal modes of operation. The most attractive 
feature of analog simulation technique is the inherent speed of 
an analog computer which makes the parametric system studies 
practical. The major drawback is in topological restructuring 
of the computer. Availability of high speed digital computers, 
the use of analog computers became obselete. The appealing 
advantage of bread board simulation turns out to be that it uses 
the system description familiar to the engineer i.e. schematic 
diagram. But this technique is rather inflexible. The parity 
simulation is also not without shortcomings. The chief among them 
are perhaps, the dedicated natural synthetic components and 
the realization are also not simpleC14ZI. 



CHAPTER 3 

MODELLING OF GTO CONVERTER 


3.1 INTRODUCTION 

Forced commutation circuitry provided with thyristors 
increases the complexity of the convereter configuration- In 
recent times, GTOs are finding application where earlier forced 
commutating thyristors were being used. This results into 
simplified converter configuration due to the absence of 
commutating elements in the power circuitry. 

In this chapter a three phase six GTOs converter system 
representation, based on graph theoretic analysis is described. The 
approach is conceptually simpler and leads to efficient 
formulation of the converter equations corresponding to all 
possible modes of operation. The state equations are derived using 
network cutset matrix which can be modified conveniently depending 
upon converter conduction patterns. The topological analysis leads 
not only to the system equations but also to the expression of 
dependent variables such as voltage across GTO valves, etc. 

The results of steady state test simulations are 
presented for both modes of the GTO converter operation 
(rectifier/ inverter) to illustrate the validity of converter 
model and the capability of computer program developed. 

3.2 CONVERTER SYSTEM 



A converter system employing six GTOs is shown in Figure 

(3.1). This includes the leakage impedance of converter transformer 

(R tL ) and the dc smoothing reactor (R.,L.). The effect of the ac 
s s d a 

system on the converter is represented by the source voltages 
* e 1 ,e 2’ e 3*‘ "^ e s y stem also includes a transient suppressor 

circuit (R t C) connected across the supply. In parallel to each GTO 
, a snubber circuit is connected comprising of a series RC 
circuit with a diode in parallel. In series with each GTO, a 
diode with R^— combination in parallel, is provided to give the 
reverse voltage blocking capability to the GTO. The dc network is 
also represented as a voltage source (V^) the converter 

terminals. 

If the objective of the digital simulation is to 
investigate the control system response in the frequency range 
upto 1 kHz, the converter circuit of Figure (3.1) can be 
simplified as follows: 

1. The snubber circuit can be neglected, since the high 
frequency switching transients are not of interest. 

Z. The diodes and R d ~C d combination in series with GTO can be 
ignored as GTO is assumed to have high reverse voltage capability. 

Based on these assumptions, the simplified converter; 
system configuration is shown in Figure (3.Z). 

3.2,1 Firing Scheme 

The output of GTO converter can be varied continuously 
using pulse width control. Pulse width is controlled between 













maximum 

and 

minimum 

values, keeping 

its 

position fixed. Among 

various 

PWM 

schemes 

the 

Equal Pulse 

Width 

Modulation (EPWM) 

is 

used 

due 

to certain 

advantages . 

EPWM 

firing scheme can 

be 


impl imented using two types of gating patterns. 

1. Symmetrical and 

2. Unsymmet ri cal 

The gating patterns are shown in Figure (3.3). In case 
of unsymmetrical gating , GTOs in upper valve group (G^G^G^) 
are gated by more than one gate pulse while those in the lower 
valve group ( Gg, , G^ , G^ ) are gated only once during a cycle. In 
symmetrical gating the switching numbers are same for all the 
GTOs. Comparing both these gating patterns, it may be observed 
from Figure (3.3) that number of switchings per ac cycle are 27 
for unsymmetrical and 24 for symmetrical gating. Therefore, 
symmetrical gating results in less switching loss which is 
uniformly distributed among all the GTOs. To avoid any switching 
transient, there should not be any gap between OFF & ON gate pulse 
applied to outgoing and incoming GTOs. 

3.2.2 Effect of Source Inductance 

In the presence of source inductance, an over-voltage 
appears across the snubber capacitance when current is interrupted 
by a GTO while turning off. The over-voltage depends on the value 
of source inductance and the magnitude of current interrupted. 
Increasing the size of snubber capacitor in order to suppress the 
over-voltage may not be appropriate as it would also increase the 



(a) UNSYMMETR1CAL GAT I N 6 



FIG.3.3 . GATING PATTERNS 













losses in the snubber circuit . A capacitor C is used across the 
source terminals to take care of the transients due to source 
inductance as shown in Figure (3.1). The value of C should not be 
very high as it would load the source permanently. A resistor R is 
provided for damping. 

The presence of source inductance also causes reduction 

in the average output voltage of the converter. Assuming that G^is 

turned off and G 3 is turned on at an instant (refer Figure 3.2). 

The current in phase B (i fa ) starts increasing from zero and would 

assume the value I . after some time. Since this current cannot 

d 

increase instanteneously and remains constant, a current i^ 

(=I d -i fa ) must be supplied by voltage absorbing circuit until i fa 

equals I d . During this process a voltage L s .(di fa )/dt appears 

across inductance L . 

s 

The average output voltage (E^) of the converter 
before the smoothing reactor can be calculated considering 
voltage drops due to the transformer resistance, inductance and 
switching devices as: 

E = E — ( 3P 6>L +XR )I -V_ (3.1) 

d O — s s d T 

Where E Q is the average output voltage when there is no drop. V T 
is drop across the switching devices under conduction and X is a 
constant. P denotes the number of pulses per ac cycle. The details 
of equation (3.1) are given in Appendix I. If R g and V ^.are assumed 
to be zero, the equation (3.1) can be written as: 



(3.E) 


E = E_-(3P)o>L I . 
d O — s d 

n 

A similar equation can be obtained for thyristor 
converters C13 ass 

E = E_— ( 3/n )c*>L I. (3.3) 

d 0 s d 

It is evident from equation (3.2) that the voltage drop 
is dependent on pulse number in addition to the transformer 
leakage reactance and dc current . Hence, the value of P shuold be 
judiciously chosen. 

3.3 MODEL FORMULATION 

The graph of Figure (3.2) is shown in Figure (3.4). This 
has 5 nodes and 12 elements. Elements 1 to 6 correspond to GTOs 
(G^ to G^). The elements 7 and 9 are equivalent representation of 
the ac system feeding the converter. Element 8 represents the dc 
network which includes dc transmission line. Elements 10 to 12 
represent the RC transient suppressor circuit. In general, the 
choice of elements considered in the formulation of a tree is 
based on the following priorities C53. 

1. Voltage Sources 

2. Capacitive Elements 

3. Resistive Elements 

4. Inductive Elements 

5. Current Sources 

For the graph shown in Figure (3.4), 4 elements 

(branches) are required to constitute a tree which for any network 





is not unique. The tree of the graph is chosen to include 


elements 10,11 and any two elements corresponding to the 
conducting GTOs , one from the upper valve group and other 
from the lower valve group. This is based on the assumption that 
the dc link current is continuous, hence one valve in the 
upper and lower valve group would always be conducting. 
Figure(3.5) shows the tree for the graph of Figure (3.4) with 
valves 6 and 1 assumed to be conducting. For any network, the 
branch and link variables are related by! 

(3.4) 

(3.5) 

where subscripts B and L of state vector refer to the branches and 
links respectively and matrix CB L 3 is the component of the 
fundamental cutset matrix. The matrices corresponding to the 
various trees which may be formed depending upon the converter 
conduction patterns, are given in the Appendix II. The branch 
variables can be further partitioned as! 


— B CB L 3 lL 
V CB i?\ 


— B = 


-B1 


V 


B2 


±B = 


-B1 


-B2 


(3.6) 


Where subscript *B1* corresponds to elements 10 & Hand *B2* 
corresponds to elements and S,, i.e. conducting GTOs included in 
the tree. Similarly, the link variables can be partitioned as 
given below: 



(3.7) 



where subscript 'L^* corresponds to elements 7,8 & 9 

*Lg* corresponds to element 12 

•Lg* corresponds to nonconducting GTOs & 

*L^* corresponds to conducting GTOs in link. 
From the partitioning of the branch and link variables, the 
matrix CB^3 can be written ass 



(3.8) 


Using equations (3.6) to (3.8), equations (3.4) & (3.5) can nov 
be written as: 


-B1 CB L11 :I -L1 CB L1Z J ±LZ CB L13' , iL3' 

CE L14 :3 -L4 

(3.9) 


CB L24 :I -L4 

(3.10) 

— L1 =BB L1 1^ T — B1 +BB L21 bT ~B2 


(3.11) 

V =CB +CB 

-L2 Le L12 J -B1 LB L22 J -B2 


(3. 12) 

v, _=cb, ^] t v„,+[:b. 

-L3 L13 — B1 L23 -B2 

- 

(3.13) 

V =EB +CB 

-L4 L14 -B1 L24 -B2 


(3.14) 


In order to simplify the model formulation, 
assumptions are made. 


the following 



1. GTOs offer zero impedance when on and infinite impedance 


when of f . 

2. Capacitances used in transient suppressor circuit are 
lossless . 

Based on the first assumption the following relations are 
evident . 

” B 2=° ' V L4 =0 ; 1 L3 =0 ».13> 

Since there is no interaction between the elements included in set 


Bg and element 12 under different operating modes of converter 
,matrix CB , ^ 3 is a null matrix. 

Thus, using equation (3.15), equations (3.9) to (3.14) 
are simplified as s 


— B1 


(3.16) 

— B2 = 

-CB LH 1 ^ L 1 ~ E " B L24 :I — L4 

(3.17) 

II 

v 

-J 

>1 


(3.18) 

” L 2 = 

“liz^bi 

(3.19) 

— L3 = 

“liz^bi 

(3.20) 

0 = 

CB L14 :l -B1 

(3.21) 


The constituent equation for elements 7,8 & 9 is given by, 

-L1 = CR a 3 ^1 +C " L a 3p -i|_1 + - (3.22) 

where 


^ L r[ V 7 ' V 8 ’ V ?] ! -LI [*7 -<8 '<•?] 


Assuming the transformer resistance and leakage reactance to be 
equal in all the three phases. 
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Where e^e^and e 3 are the three phase ac 
voltage source representing dc network, 
resistance and inductance respectively, 
basic equation is! 


P v B1 = CR f 3pi B1+ ccpi B1 


o = 2nf 

supply voltages, is 

and are smoothing 

For branches lO & 11, the 

C3.24) 


where 



I 

o 

£ 


' 1 0 ' 

CR f 3= 


? ccp= 

oC 

A 


O R 


0 

6)C 


;R and C are transient 
suppressor circuit 
elements . 
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' i io ' 


« i ss 

’ — B1 


v n 


*11 


Similarly for element 12, the following expression can be written. 

P -L2 =R ' P -L2 + C fi|_2 * C f «1/«C (3.25) 

With suitable manipulations of equations (3.16) to (3.25), the 

following equation is obtained: 




( 3 . 26 ) 


The detailed derivation of equation (3.26) is give in 
Appendix III. The formulation of state equation for two different 
modes of converter operation is illustrated using the generalised 
state equation (3.26). 

CASE I 

This case corresponds to the power mode operation of the 
converter in which say GTOs G^ and G^ are conducting. The 
submatrices for this mode of operation can be derived based on the 
cutset matrix formed for the tree with 5^=6 and S 2 =1, as given 
below: 





CK,_ 3= — ( 1 / ( 3R ) ) « CB. „„3 T CR-3CB, . .3 
5 L12 f L11 


=^-(1/3) (1/3) (1/3) J 

CK 6 1=0 ; CKg3=0 

CK 7 3=(1/(3^RC) ) . £-1 1 1 J 

CK 11 3=C03 ; LK ^2=1101 ; 0(^3=1103 ; 0(^3=1103 


CK 19 3-r[|3 l11 3*CB l1e 3CK 3 3] 



(1/3) (2/3) (2/3) 

(2/3) -(1/3) -(1/3) 

LK nn 3=L03 ; CK 22 3=C:03 ; LK^1=L03 


'20 


16' 


C K 2i 3 


(1/3) (2/3) (2/3) 

-(2/3) -(1/3) -(1/3) 


63 C 


CK 15 3=CK 17 aCK 1 D+CK E0 aCK 2 3+CK E1 3 


C k 17 3= 


-1 


L 6>L 


0 
1 

iiiL J 
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(R/3) 
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1 tiiL i^RCj 
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uL 


-2 

uL 


6lL 
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(R/3) 


s J d 

The final set of equations can be obtained by putting above 
calculated submatrices in equation (3.26). The submatrices of 
different GTO conduction patterns can be calculated in the same 
manner. 


CASE II 



This case corresponds to the free wheel mode of converter 


operation in which S^= 1 and Sg= 4. The submatrices for this mode 
of operation can be calculated as followss 


The value of submatrices CK.I1 and CK_I1 will be same as in 

1 3 


case I, 



' 0 

0 

1' 


-I" 


' 1 

-1 -1 

1" 

CB. ..3= 
L1 1 

-1 

0 

0 

; CB, „□= 

* L12 

-1 

’ ^LIS 3 

1 

0 

1 

-A. 

0 


“LEI 3 " 


O -1 o' 
0-1 0 


‘ CB LE3 :1 - 


0 110' 
10 0 1 


? CB L14 D=CB L24 3=n03 


CK_1=C6jL □ ^CB. j|j , 3 T 
2 a L1 1 


‘2 0 -1 ' 


' 0 -1' 


-1 -2 

3ooL 3oiL 




3ti>L 3a)L 

s s 




s s 

o 

T* , 

o 


o 

o 

— 

0 0 

“ L d 





- 1 2 


A A 


2 1 

3toL 3uL 

s s J 


A U 


363 L 363L 

s s J 


CK^U=— ( 1 / ( 3R ) ) .CB L12 3 CR f 3CB L11 3 


=£-(1/3) 0 (1/3)J 

CK 7 3=[-1 O 1 ] .^ c 


CK,3=C03; CK_3=i:03j CK . „ 3=C0:i ,* 

o o n 

CK 12 D=C03; CK 13 3=tI03; CK 14 M03 


CK 2Q 3=C03; CK 22 3=C03 ; HK^^COU 



10 2 


. (R/3) 

19 -2 0 -1 



CK 15 3= CK 19 3CK i n+CK 20 3IIK 2 3+CK 21 3 


f -R 

, 1 1 


f — 2R 

2 1 

s 

61RC J 

0 

[ uL l 

+ <*RC J 


(R/3) 


f 

- 2 1 

A 

r ^ 

1 1 

l “L l 

oRC J 

U 

[ 6JL 

6>rc J 


1 0 * 

,(R/(3wL s ) ) 

0 1 

0 0 1 

10 0 

The final equation can be obtained by putting these calculated 
values of submatrices. For other GTQ conduction patterns the value 
of submatrices can be calculated in the same way. 

DISCUSSION 





The equation (3.26) is formed at each time step. The 



changes in the status of GTOs in the link can be considered by 
rearranging the columns of matrix CB L 1. The matrix HB^3 has to be 
changed as the status of values in branch changes. The method given 
here is simpler, both conceptually and computationally. Also the 
GTO characteristics can be assumed different than that of an 
ideal switch. The constituent equation in this case would be 
different than equation (3.15). To give an illustration, the 
impedence of the GTO during nonconduction can be assumed to be 
finite (determined by the grading and damping circuit). Inclusion 
of this characteristic can be easily accomplished. In contrast, 
the ref erenceC153 implicitly assumes the GTO valve impedence to be 
infinite during nonconducting period. The snubber circuit may 
also be included but then complexity of graph will be increased. 
Inclusion of snubber circuit will increase nodes and elements of 
the graph. In contrast, the reference C3U has implicitly neglected 
the snubber circuit. The number of state equations per converter 
vary depending upon the number of conducting GTOs (N) and is given 
by (N+g). 

3.4 COMPUTER PROGRAM 

Based on the procedure outlined in the previous : 
sections, a computer program has been developed to validate the j 
converter model. The structure of the program is modular with 
subsystem described in the individual subroutine. At a particular 
instant of time, the states of converter are defined and the 
equations are formulated using the converter model described 



earlier. The converter state changes due to the extinction of the 
GTOs or the firing of the GTOs. The converter state is checked at 
the begining of each integration time step. The presence of firing 
pulse indicates the conducting GTO which is generated by pulse 
generating subroutine and the GTO is put into the conducting state 
if it is high and the GTO is put into nonconducting state if it is 
zero. The salient feature of the program is given below. The flow 
chart of the program is given in Figure (3.6). 

INITIAL CONDITIONS 



The simulation 

can 

proceed from a) 

zero 

initial 

conditions or b) 

initial conditions established 

from the 

steady 

state 

operating 

point. For 

the 

second case the 

requi red 

steady 

state 

operating 

conditions 

are 

calculated. The 

present 

program 


uses the EPWM firing scheme. The program can be run more than a 
cycle. The voltage of ac supply, frequency of supply, transformer 
leakage impedance , smoothing impedance, step of integration, 
output voltage pulse per ac cycle, modulation index and reduced 
bus incidence matrix have to be given as the input datas. The 
analysis is started assuming that the GTOs G^ and Gg are in 
conduction prior to turning on of G^ at a which is calculated by 
expression given in table (2.1) for EPWM. The conduction pattern 
of the GTOs are given in the Appendix V for symmetrical gating of 
EPWM scheme. The initial conditions required for the numerical 
integration can be fed externally or determined by the program 
using following relations. 















Initial load current as 


I d ( 0 )=i 8 (°)= o.a.(E d -E b >/R i 

inhere E Q is the average output voltage and can be calculated by 
lsing equation (1.4). Other initial conditions can be calculated 
as follows. 


i 7 (0)= 0.0 
i 9 <0)= 0.0 
V 10 <°>= e f e 


V^CO)- e 2 -e 
3.5 TEST SIMULATIONS 


2 

3 


The capability of the program developed is illustrated 
by simulating the GTO converter with an active load. The details 
of the system parameters are given in Appendix IV. The test 
simulations are carried out to investigate the system response in 
steady state for both modes of converter 

( rectification/ inversion) . 

Figure (3.7) shows the steady state output current of 
rectifier and the ac phase *A* current is shown in Figure (3.8). 
The valve voltages across G^ and G,, are given in Figure (3.9) to 
Figure (3.10). Figures (3.11) to (3.13) show the current through 
GTO valve 1, current through transient suppressor circuit and dc 
output voltage of rectifier operation. The voltage across transient 
suppressor circuit elements are given in Figures (3.14) & (3.15). 

Figures (3.16) to (3.21) show the result obtained from 
inverter operation of GTO converter. The input dc voltage applied 
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FIG. 3.15 VOLTAGE ACROSS TRANSIENT SUPPRESSOR CIRCUIT 



to the inverter is 180.0 volts with negative polarity and other 
parameters remain same. Figure (3.16) shows the ac output current 
from inverter. The waveform of dc current flowing through this is 
shown in Figure (3. 17). The current through GTO valve 1 and 
suppressor circuit is shown in Figures (3.18) and (3.19) 
respectively. The voltage across GTO valve 6 and across transient 
suppressor circuit are also given in Figures (3. HO) and (3. HI) 
respectively. 

The harmonic analysis of rectifier operation of GTO 
converter is also presented under assumption mentioned in previous 
sections. This is done by discrete fourier analysis. The 
X-Axis of graph represents the multiple of fundamental frequency. 
The Y— Axis shows the magnitude of harmonics. Figure (3.HH) shows 
the harmonic components of ac voltage and Figure (3.H3) shows the 
harmonics in ac phase current. The harmonics components of dc 
sides voltage and current are shown in Figures (3.24) and (3.25) 
respectively. 

DISCUSSION OF RESULTS 

The results demonstrate the operation of converter 
(rectifier / inverter) under steady state. In Figures (3.7) and 
(3.16), ^ ac current of 1 ampere (approx.) flows continuously 
which leads the phase current by an angle of 90°. This is only due 
to the RC transient suppressor circuit, present at converter ac 
bus. The value of C must be chosen judiciously so that it should 
not load the ac source substantially. The value of R and C can be 
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FIG. 3.23 HARMONIC COMPONENTS OF RECTIFIER AC CURRENT 







optimized accordingly but this optimisation is not attempted here. 
The spikes present are due to change over of GTO conduction 
pattern which causes the interruption of source current. The dc 
current has ripples. The current rises in power mode and falls in 
free wheeling mode. This rise and fall will depend on the value of 
resistances and inductances present in the dc current path. Apart 
from fundamental component in ac voltage are not of significant 
magnitude. The hormonics present in source current are 2nd to 12th 
however, the magnitudes of 5th ,7th, 10th and 12th are higher than 
others. The magnitudes of harmonics are dependent ofs 

1. method of modulation and number of pulses per cycle of 

* 

source voltage. 

2. nature of loads 

3- presence of source impedance and 

4. presence of RC transient suppressor circuit. 

The harmonic in dc side is only 11th, although 10th, 12th and 13th 
are also present but their magnitudes are lower than 11th. 


3.6 CONCLUSION 


In this chapter a state space model for the three phase 
ac-dc PWM GTO converter system has been developed based on the 
graph theoretic approach. The generalise set of state equations 
has been derived based on the network cutest matrix. The system 
model is general enough to be used for investigation of 
performance for three or more GTO conduction at a time. A 
simulation program bas^ybeejn , wr-H^^F\ to validate the model 


Ace. jNo. 





developed. The 
parameters for 


test simulation has been carried out 
12 pulse per cycle of source voltage. 


with 


given 



CHAPTER 4- 

GTO INVERTER IN HVDC LINK 


4.1 INTRODUCTION 

An HVDC system, with thyristor rectifier and their 
associated controls, dc network consisting of transmission line 
and GTO inverter with controls are discussed in this chapter. 
Digital simulation requires adequate representation of the various 
components. This chapter describes the modelling of the 
constituents of an HVDC system in brief. Emphasis is laid on 
interfacing of GTO inverter and its controls in place of thyristor 
inverter using appropriate interfacing variables. A suitable model 
for GTO has been developed in previous chapter. These models are 
used for digital simulation to study the performance of HVDC link 
with GTO inverter. The simulation has been carried out with ac 
system representation. A computer program is used to carry out a 
large number of test simulations, both under steady state and 
transient conditions, to demonstrate the feasibility of GTO 
inverter in two terminal HVDC link. The disturbances considered, to 
study the dynamic behaviour, are change in current reference 
setting and ac voltage dips at both rectifier and inverter 
terminals. 

4.2 CONVERTER REPRESENTATION 


A six pulse thyristor 


converter is modelled as 


a 



variable voltage source behind a variable impedance. Reference 

C173 describes the development of an equivalent circuit for a six 

pulse thyristor converter. This is based on graph theoratic 

approach which enables to formulate the converter equations for 

all possible modes of operation assuming continuous dc link 

current. Converter equations are derived using cutset matrix and 

the basic circuit equations. The voltage source (e ) in the 

eq 

equivalent circuit is a function of ac bus voltage and hence has 

to be computed at each time instant. The equivalent circuit 

parameters R and L are dependent on the converter conduction 

pattern and are recalculated every time when the conduction 

pattern is changed. The V c has been included to represent the 

effect of the dc system. Both ac and dc voltage sources are not 

constant. R . and L. denote the resistance and inductance of 
d d 

smoothing reactor. The thyristor converter equivalent circuit is 
shown in Figure (4.1). The thyristor converter is used as 
controlled rectifier. 

The GTO converter is used as inverter in HVDC link. The 
detailed representation of the GTO converter is given in chapter 
3. 

4-3 DC SYSTEM REPRESENTATION 

The dc system comprises of the filters and transmission 
line. The transmission line is modelled as n- equivalent circuit 
with the shunt arms consisting of capacitors and the series arm a 
combination of resistance and inductance. A more practical 




FIG. 4.1 EQUIVALENT CIRCUIT OF 6 PULSE 
CONVERTER 



representation of the line is a number of n— sections. This is 
achieved by dividing the line into a number of sections, 
representing each section by its own n— network and connecting all 
such networks in series. In the dc system representation, the 
converters are represented as current sources as shown in Figure 
(4.2). The dc filter used for filtering any harmonics present in 
the dc current at rectifier end is represented in the dc system 
model. The dc filter at inverter end is not represented. 

4.4 RECTIFIER END CONTROLS 

Generally, rectifier is equipped with a constant current 
controller. The control of the rectifier considered here is based 
on digital control technique proposed by Freris et al C161. The 
firing scheme is basically an Equidistant Pulse Control (EPC) with 
pulse frequency control (PFC). The interfiring period (IFP), which 
is the interval between two successive firing instant, is 
calculated as 

IFP « 60°+ G i 

Where G^ is the firing correction, calculated by controller . In 

steady state, Gj=0 and firing takes place at 60° intervals for a 

six pulse rectifier. Figure (4.3) shows the block diagram for 

calculation of G. and IFP for constant current control. In the 

1 

case of constant current control 

Q.= K.V 
1 c 

where K is the gain and V c is the control signal. The control 
signal V c may be a linear function of current error, in which case 







Arrillaga C17H, however, has suggested that the 
inclusion of a derivative term in the calculation of control 
voltage improves the dynamic behaviour of the system. Then the 
control signal V c becomes 

V c= K 1* (I d- I ref )+ K 2 (dI d /dt) 

This control signal V £ is sampled at an appropriate time and this 
sampled value forms the basis for the calculation of Qj . 



FIG. 4-3 CONVERTER CONSTANT CURRENT CONTROL 


The other control at rectifier end is minimum firing 
angle control. In this control, the implementation technique is to 
calculate the voltage across the valve corresponding to the 
minimum firing angle and if at the instant of firing, the voltage 
across valve is less than this voltage, the valve conduction is 
prevented till the voltage rises to the rquired level. The 
transition between minimun firing angle control and constant 
current control is thus automatic C203. 





4.5 INVERTER END CONTROLS 


The thyristor inverter control was constant extinction 
angle control and constant current control. Under normal 
conditions, the operation was in constant extinction angle control 
mode. Since GTO thyristor is used at inverter end and there is no 
cocept like extinction angle because of the facility of extinction 
of GTO valve at any desired instant. The firing time and 
extinction time of a GTO converter is decided by modulation index 
which can be varried between O and 1. Since regeneration time is 
required for normal operation so the value of modulation index is 
limited between 0.05 to 0.95. In this thesis, EPWM scheme with 12 
voltage pulses per ac cycle is used. The firing and extinction 
instants of GTOs are calculated by expression given in table 
( 2 . 1 ) . 



FIG. 4.4- BLOCK DIAGRAM OF INVERTER CURRENT CONTROL 
Basically there are three terms in the expression (1.4) 
of Appendix I (i.e. dc voltage, dc current & modulation index). 





Therefore one can control any of them or two of them 

simultaneously or even all of them. Since conventional thyristor 

inverters have constant current control and constant extinction 
angle control (CEA) . In the same way here performance study is 
tried for GTO inverter using one of the following two 
alternatives . They are : 

d-Jt s/CitTr 

1. constant current control and constant angdse 

control 

2. constant current control and constant modulation index 
control . 

In PWM GTO converter, two modes of operation is 

performed, one is power mode in which ac source is directly 

connected to the dc system and other is freewheeling mode in which 
ac source is disconnected and there is dc short circuit. Due to 
large variation in dc voltage, the use of constant voltage 
control is very difficult, however, voltage transducer may be used 
to solve this problem but it will increase the system time 
constant. In this thesis, the study of the GTO inverter control is 
performed on the constant current control and constant modulation 
index. Constant modulation control is similar to the CEA control 
of thyristor inverter. In order to avoid the operation of both 
rectifier end current controller and inverter end current 
controller, a margin of 1051 in current refference is given. When 
the inverter dc current decreases below 0.9 p. u . , the current 
controller of inverter end starts working and for other cases 



modulation index control. The value of margin is so chosen that 
the transition beween inverter end controller is smooth. The 
output of constant current controller is directly related to the 
modulation index and firing & extinction angles are calculated at 
every instants. The block diagram of constant current control of 
inverter has been given in Figure (4.4). 

4.6 COMPUTER PROGRAM 

The development of a digital computer program for 
simulating an HVDC system with thyristor rectifier and inverter is 
described in reference C171. This incorporates both 6 and 12 pulse 
converter representation along with firing control schemes based 
on IPC and EPC. The program is extremely modular and therefore 
permits implementation of different control schemes and other 
advance features necessary to simulate a practical system. The 
program is augmented by Ajai C20I1 to include digital control 
scheme, dc filters and dc line faults. This program is further 
augmented to replace thyristor inverter by GTO inverter. Suitable 
control schemes are incorporated for the simulation of an HVDC 
system with GTO inverter. The basic flow chart is same as in 
referenceC203. The computation of modulation index from controller 
output, hence firing and extinction instant of GTO valves are 
calculated to check the incoming and outgoing GTOs. The program 
which was augmented and reported in reference C20U is used 
suitably with some modification for GTO inverter in place of 
thyristor inverter. 



The simulation can start either assuming the system 
under steady state operating condition or with initial condition. 


In the case of 

f ormer 

the initial 

conditions 

for 

various state 

variables 

have 

to 

be supplied. 

The calculation 

of 

initial 

conditions 

for 

all 

state variables except 

for 

GTO 

inverter 

variables 

have 

been 

reported in 

reference 

C203 . 

The 

initial 


conditions for GTO inverter state variables are calculated in 
chapter 3. 

4.7 CASE STUDY 

The effect of GTO inverter in HVDC system is 
investigated in detail by digital simulation. The operating 
conditions and parameters are given in Appendix VI which is taken 
from reference CHOU. The control parameters for GTO inverter are 
also included. The detail investigation includes the comparision 
of rectifier / inverter dc voltage and current with thyristor 
rectifier / inverter , harmonic analysis and reactive power 
requirement. 

4.7.1 Steady State Simulation 

To demonstrate the program capability, steady state 
performance of HVDC system with six pulse converter conf iguration 
at rectifier end and GTO converter at inverter end has been 
studied. Figures (4.5) to (4.8) are the instanteneous plots of 
rectifier dc voltage, dc current, inverter dc voltage and dc 
current respectively. The average value of rectifier dc voltage 
and dc current, inverter dc voltage and dc current are given in 
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FIG. 4.5 RECTIFIER INSTANTENEOUS DC VOLTAGE UNDER STEADY 
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FIG. 4.7 INVERTER INSTANTENEOUS DC VOLTAGE UNDER STEADY STATE 
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FIG. 4.8 INVERTER INSTANTENEOUS DC CURRENT UNDER STEADY ST, 



Figures (4.9)to (4. IE). The harmonic components of ac and dc sides 
of inverter end are shown in Figures (4.13) to (4.16). The ac 
current in a phase and ac voltage are given in Figure (4.17) to 
show the power factor. The harmonics and power factor of thyristor 
inverter are well known C13. 

4 .7.2 Transient Response 

The dynamical behaviour of the HVDC system employing 6 
pulse thyristor rectifier and GTO inverter has been investigated 
with current reference change and ac voltage dip on the both sides 
of converter. The average dc voltage and average dc current at 
both rectifier and inverter end, harmonic components of ac and dc 

values (voltage/current) of inverter end during fault( 1st cycle) 

/ 

are plotted. The following cases are studied. 


Case No. Name of case 


1. Change in current reference by 20% 

2. 99% dip in one phase for 10 cycle at rectifier end 

3 . 99% dip in two phases for 5 cycle at rectifier end 

4. 50% dip in all phases for 5 cycle at rectifier end 

5 . 99% dip in one phase for 10 cycle at inverter end 
99 % dip in two phases for 5 cycle at inverter end 

7 . 99X dip in all phases for 5 cycle at inverter end 

8. 50% dip in all phases for 5 cycle at inverter end 



Idcr 



FIG. 4.10 RECTIFIER AVERAGE DC CURRENT UNDER STEADY STATE 
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FIG. 4.13 HARMONIC COMPONENTS OF INVERTER AC VOLTAGE 
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FIG. 4 . 17 (a) INVERTER PHASE X' CURRENT 
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(1) Change in Current Reference: 

The response for two terminal case shows that the 
current controller at rectifier end comes into action as a step 
change in the current refernce is applied to the both ends. The 
inverter current controller will work only when dc current will 
fall down below set the value. The application of reference change, 
the current settles down to the new value in certain time which is 
dependent on gain and time constant of controllers. Figure (4.18) 
shows the rectifier end average voltage and current 

following 2 . 07 . increase in current reference at both ends of the 
converter. The inverter end average dc voltage and dc current are 
shown in Figures (4.19) 4 (4.20) respectively. 

(2) . 997 . dip in one phase at rectifier: 

The probability of this fault is very high but severety 
is low. Figures (4.21) 4(4.22) show the average dc voltage and dc 
current at rectifier end, and inverter end average dc voltage and 
current are given in Figures (4.23) 4 (4.24) respectively. The 

harmonic components of inverter end ac voltage, ac current, dc 
voltage, dc current are given in Figures (4.25) to (4.28) 
respectively. The power factor can be seen by Figures (4.29) and 
(4.30). 

(3) 997. dip in two phases at rectifier : 

The dip of 99 X is provided in phases 'A’ 4 *B* for 5 

cycle at the begining of 5th cycle. The average dc voltage and 
current of rectifier terminal are shown in Figures (4.31) 4 
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FIG. 4..18(b) RECTIFIER TERMINAL DC CURRENT FOR CASE 1 
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FIG. 4.21 RECTIFIER TERMINAL DC VOLTAGE FOR CASE 2 



FIG. 4-.es RECTIFIER TERMINAL DC CURRENT FOR CASE 2 
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FIG. 4,25 HARMONIC COMPONENTS OF INVERTER AC VOLTAGE FOR CASE 2 
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PIG. 4£6 HARMONIC COMPONENTS OF INVERTER AC CURRENT FOR CASE 
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FIG. 4.27 HARMONIC COMPONENTS OF INVERTER DC VOLTAGE FOR CASE 
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FIG. 428 HARMONIC COMPONENTS OF INVERTER DC CURRENT FOR CASE 
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jrjQ 4. 32 RECTIFIER TERMINAL DC CURRENT FOR CASE 3 


(4.32) « Similarly, the inverter dc voltage and current are shown 
in Figures (4.33) & (4.34). The harmonic components of inverter ac 
voltage, ac current, dc voltage, dc current are shown in Figures 
(4.35) to (4.38) respectively. 

(4) 507. dip in all phases at rectifier 

The Figures (4.39) to (4.42) show the rectifier dc 
voltage & dc current , inverter dc voltage & dc current 
respectively. The harmonic analysis of inverter ac voltage, ac 
current, dc voltage, dc current are given in Figures (4.43) to 
(4.46) . 

(5) 99 7. dip in one phase at inverter end 

Figures (4.47) and (4.48) show the rectifier average dc 
voltage and dc current. The invreter average dc voltage and dc 
current are shown in Figures (4.49) and (4.50). The harmonic 
spectrum of inverter ac voltage, ac current, dc volatge, dc 
current are given in Figures (4.51) to (4.54). Due to the fault 
2nd and 3rd harmonic components are generated in the dc current 
and 2nd harmonic in dc voltage. 

(6) 99% dip in two phases at invereter 

The waveforms of average dc voltage, dc current at 
rectifier terminal and average dc voltage, dc current at inverter 
terminal are shown in Figures (4.55) to (4.58) respectively. The 
harmonic spectrum of inverter ac volatge, ac current, dc voltage, 
dc current are given in Figures (4.59) to (4.62). During fault the 
2nd and 3rd harmonics are generated in the inverter ac voltage and 




FIG. 4.33 INVERTER TERMINAL DC VOLTAGE FOR CASE 3 
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FIG. 4.35 HARMONIC COMPONENTS OF INVERTER AC VOLTAGE FOR CASE 3 
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FIG 4.36 HARMONIC COMPONENTS OF INVERTER AC CURRENT FOR CASE 



MAGNITUDE 




o.co 

- 1.00 1,00 3,00 5,00 7.00 S. 00 11,00 13.00 15,00 17.00 19.00 21.00 

FREQ. 

FIG. 4-37 HARMONIC COMPONENTS OF INVERTER DC VOLTAGE FOR CASE 3 
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FIG. 4.38 HARMONIC COMPONENTS OF INVERTER DC CURRENT FOR CASE Z 
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FIG. 439 RECTIFIER TERMINAL DC VOLTAGE FOR CASE 4 
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FIG. 4.40 RECTIFIER TERMINAL DC CURRENT FOR CASE 4 
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FIG. 4.43 HARMONIC COMPONENTS OF INVERTER AC VOLTAGE FOR CASE 4 
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FIG. 4.45 HARMONIC COMPONENTS OF INVERTER DC VOLTAGE FOR CASE 4 
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FIG A.AR HARMONIC COMPONENTS OF INVERTER DC CURRENT FOR CASE 4 






FIG. 4.47 RECTIFIER TERMINAL DC VOLTAGE FOR CASE 5 



FIG. 4.48 RECTIFIER TERMINAL DC CURRENT FOR CASE 5 



FIG. 4.4-9 INVERTER TERMINAL DC VOLTAGE FOR CASE 5 
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FIG. 4.50 INVERTER TERMINAL DC CURRENT FOR CASE 5 



FIG. 4.51 HARMONIC COMPONENTS OF INVERTER AC VOLTAGE FOR CASE 5 
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no. «56 RECTIFIER TERMINAL DC CURRENT FOR CASE 6 



FIG. 4.57 INVERTER TERMINAL DC VOLTAGE FOR CASE 6 



FIG. 4J58 INVERTER TERMINAL DC CURRENT FOR CASE 6 
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FIG. 4.59 HARMONIC COMPONENTS OF INVERTER AC VOLTAGE FOR CASE 6 



FIG. 4.60 HARMONIC COMPONENTS OF INVERTER AC CURRENT FOR CASE 6 
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FIG. 4,61 HARMONIC COMPONENTS OF INVERTER DC VOLTAGE FOR CASE 6 
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FIG. 4.62 HARMONIC COMPONENTS OF INVERTER DC CURRENT FOR CASE 6 




ac current. 

(7) 99 7. dip in all phases at inverter 

Figures (4.63) to (4.66) show the rectifier average dc 
current, dc current and inverter average dc voltage, dc current 
respectively. The harmonic components of inverter ac voltage, ac 
current, dc voltage, dc current are shown in Figures (4.67) to 
(4.70). This fault is very severe but less probable. This is 
studied for 5 cycle dip only. 

(8) 50% dip in all phases at inverter 

This fault is known as three phase remote fault. This 
effect is rather less at inverter than 997. dip at inverter. 
Figures (4.71) to (4.72) show the rectifier average dc voltage, dc 
current and inverter average dc voltage, dc current respectively. 
The harmonic components of inverter ac voltage, ac current, dc 
voltage, dc current are shown in Figures (4.73a) to (4.74b). 
DISCUSSION OF RESULTS 

The steady state behaviour of HVDC system with GTO 
inverter is improved than using GTO converter with an active load. 
The lower order harmonics are totally eliminated. The dc output 
voltage have only 11th harmonic component. During disturbances, 
particularly dips in ac voltage at rectifier end, no extra 
hormonic components are generated in inverter side quantities. 
Some harmonics are generated during faults at inverter ac voltage. 
The power factor is maintained at unity. The comparision of HVDC 
system with GTO inverter has been made with thyristor inverter. 





FIG. 4- 65 INVERTER TERMINAL DC VOLTAGE FOR CASE 7 



FIG. 4-.6S INVERTER TERMINAL DC CURRENT FOR CASE 7 
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FIG. 4.68 HARMONIC COMPONENTS OF INVERTER AC CURRENT FOR CASE 7 
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FIG. 4.71(b) RECTIFIER TERMINAL DC CURRENT FOR CASE 8 
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FIG. 4.73 <a> HARMONIC COMPONENTS OF INVERTER AC VOLTAGE F< 
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FIG. 4,74 (a) HARMONIC COMPS. OF INVERTER DC VOLTAGE FOR CASE 8 



FIG. 4.74 <b) HARMONIC COMPS. OF INVERTER DC CURRENT FOR CASE 8 
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4.77 50/., THREE PHASE, 10 CYCLE DIP AT INVERTER 

NORMAL RECOVERY 




Figures (4.75) to (4.77) are the system performance with thyristor 
inverter for case 5,7,8 respectively. 

4.8 CONCLUSION 

This chapter has been devoted to study the two terminal 
HVDC system performance with GTO inverter. A control strategy is 
made for GTO inverter to co-ordinate the conventional rectifier 
control. The optimization of control parameters has not been 
attempted. Both steady state and transient behaviour of the system 
without ac system representation have been studied. A 
comparision of the HVDC system with GTO inverter has been made 
with thyristor inverter. 



CHAPTER 5 

CONCLUSIONS AND FUTURE SCOPE OF WORK 


This thesis has been mainly devoted to the development 
of a suitable model of the converter, consisting of six GTOs, for 
dynamic digital simulation. This model has been developed based on 
the graph theoretic approach which is conceptually simpler and 
leads to efficient formulation of the converter state equations 
corresponding to all possible modes of operation. The state 
equations are derived using network cutset matrix which can be 
modified conveniently depending upon converter conduction 
patterns. The investigation of performance of the GTO converter 
with active load has been made in both the modes of converter 
operation ( rectification/inversion) . 

A major area of emphasis has been on the study of the 
performance of GTO inverter in HVDC link with appropriate inverter 
controls. The following conclusions have been drawn from above 
studies. 

1. The lower order harmonics, which were present in the 
steady state operation of the converter with active load, are 
eliminated when GTO converter is used as inverter in HVDC link. 

2. There is no requirement of reactive power at inverter end 
< i .e. power factor is maintained at unity). 

3. The magnitudes of over-voltages due to fault are reduced 


14 



significantly with GTO inverter. 

4. System recovers the steady state after clearing the fault 
without much disturbances within reasonable time. 

5. The system performance is less affected to the rectifier 
faults than the inverter faults. 

Since the major emphasis of this thesis is software 
developement to see the feasibility of GTO inverter in the HVDC 
system, the optimization of the controller parameters and the, 
transient suppressor circuit parameters may lead to a further 
improvement in system performance. The following further additions 
can be made s 

1. A suitable model can be developed using snubber circuit 
and also assuming that the switching devices are not ideal. 

2. Other PWM schemes can be used to investigate the 
feasibility of GTO inverter in HVDC link. 

3. The detailed representation of ac system can be included 
at both sides of the converter, so that performance of GTO 
inverter with weak ac system can be studied. 

4. It may be thought to make leading power factor instead, of 
unity power factor by suitable PWM schemes. 
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APPENDIX I 


Voltage drop occurs only P times for a half cycle of the 

current i consists of P pulses 

c _ . l 1 . (di./dt) .dt.2.P.3 

E x “ s b 

X O s 1/f 

= S d L .di..6.P.f 
0 s b 

= L .1 ,.6.P.f 
s d 

=(3P/n) .2.rcf .L *I_, 
s d 

E = (3P/n ) ,6>L . I (1.1) 

X S u 

Where t^is a time required for the current i^ to reach the value 
I d . Average output voltage of converter, 

E d = E Q -E x (X.2) 

where E Q is avearge output voltage when L g = 0.0 
E q can be calculated as follows : 

P n/b+fi. 

E == (1/(n/3))E S 1 d(^t) (1.3) 

° i=1 n/6+«. AB 
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v P h 
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v P h 

Sin 
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E « 
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V2 

v P h 
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E 

AB 

*y6. 
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.Sin 

(iot+n/6) 



* p n / 6+/3 . 

=-4=*J6. V . .r / 1 Sin (<ot+© ) d(wt) 

n/3 ph . ... r 

r i=1 n/6+«. 

l 


< 1 .4) 


where p=P/6 and & =n/6 

r 

Solving and simplifying the equation (1.4), the following equation 


can be obtained : 

P 

E_ = V . r C Cos (nr/3+« . )— Cos ln/3+ft. >3 (1-5) 

O do i i 


where V . = 3/3 V2 V . 

do ph 


Voltage drop due to ac line resistance, inductance and switching 
devices is given below. 

E x = (3P 6 jL_ s + X.R s ).I d -V T (1.6) 

n 

where L and R are per phase source inductance and resistance 

s s 

respectively. 


Vy= forward voltage of diodes*2 + 2* forward voltage of GTO 

. __ Total conduction period of succesive ot and ft. 

A* ** * 


T= 1/(2pf) 


Hence average output voltage 



APPENDIX II 


The reduced incidence matrix of the graph given in 
chapter 3, taking node one as reference node is; 



The cut set matrices can be calculated by reduced incidence 

matrix EA3 C2-1H. 

There are two mode of operations, one is power mode and other 
is freewheeling mode. Each mode have six cutset 
matrices. The twelve cutset matrices each corresponding to 

the valve considered in the tree ,are given below. The superscript 
ij denotes the valve number , and taken in branch. 



( 

f 










6 


0-10 0O101 


For free-wheeling mode. 



□ 




□ 





same 


All the other node ere CB L 3 , CE> L 1 , and CB L 3. They 

63 14 25 

as matrices CB L 3,CB L 3 and CB L □ respectively 


are 

Only 


def f errence is interchange of last two rows 



APPENDIX III 


Putting the value of V, rom eouatjon (3.18) into 

— l_ 1 

equation (3. 14) , f ol lowing equation can be obtained : 

= [cR a 3+uCU a 3pJ i L1+ e 
Above equation can be rearranged as 

" 1l 1= - C “ L a ;3 ‘ 1 iL1 +I: “ L a^ 1,:B L 1 i :,T yB 1 -' : ‘' 

or, 

p Al1 = + CK 2- 1 -B1 + LK 3 2 - (III.1) 

where CK„3= -QiL 3 _1 CR 3 
1 a a 

CK_3= CoiL l“ 1 CB 1 _3 T 

H a L1 1 

CK„3= -CoL 3 -1 
u a 

Premultiply CB^.^3 7 in equation (3.24) and equating with equation 

(3.21) , 

CB. a „3 7 pV„ „ = CB. J .3 7 CR„3pi„ a + CB. .3 7 CCl3i_ 

L14 — B1 L14 l V f -B 1 L14 f-BI 

or, 

C03= CB L14 3 T CR f ]pi B1+ CB L14 ] T CCp iB1 
Putting the value of i^ in above equation from equation (3.16), 

CM- CB L14 3 T CR f 3p[-CB L11 ^ L1 -CB L1z 3i L2 -CB^Ili^] - 


CB 

Premultiply CB L 12 


1 

L14 

3 7 in 


T 


CCj3^-CB L11 3i L1 
equation (3.24) 


-“LIZ^LE - CB L14 ] iL4] ,II1 - Z1 
and equating with equation 


(3.19), this can be obtained 



CB L12 :i P -B1 CB L12 3 *” R f ^ P — B1 + CB L12 3 CC f 3 -B1 

and 

P[ ' B L12 :I -B1 P — L2 

- p — LE = “lIB^VIbI* CB L12 3 T CCp iB1 
Comparing above equation to equation (3.25) 

R ' P i.L2 +C f L2” CB L12' 3 C ' R f ' lp iB1 +I ' B L12 :i CC f 3 i.B1 

~ BB L1 1 11 — L1~ BB L12^— L2 +BB L14^— L 4 J 

Rearranging above equation, the following equation ca be ablained: 
[R + CB L12 3 T 2R f][ ; B L12 3] pi L2 = -[c f+ CB L12] T CCpCB L12] ]i L2 - 

The above equation can be simplified as, 

pi, = -CK.Ii, _+CK,_3pi ^CR^pi, A +£K^li, j| +CK_li, . (III. 3) 

— L2 4 — L2 5 —LI 6 — L4 7 —LI 8 — L4 

where , 

CK^3= -C^/R (always) 

CK 5 3= -(1/3R).CB L12 3 T CR f TCB L11 3 
CK 6 3= -(1/3R).CB L12 D T CR f ICB L14 I 
CK y 3= -(1/3R).CB L12 T T CC'ICB L11 3 
LK s 3= -(1/3R).CB l 12 3 T CCJ.3CB l14 3 
Equation (III. 2) can be further simplified as follows: 



‘^W^LII 31 ’ pi^ZB^ZB^JZCpi^ZB^ZRpZB^ 

PiLE +CB L14 :,T[:B L1 2 :3C fiL 2 +I:B L14 ;!T|:B L 1 4 :,i;R f :) PiL4 + 

CB. . .3 T CCI3CB, . .Hi. . 

L14 f L14 — L4 

If the values of resistances and capacitances of transient 
suppresor circuit of each phase are equal, then the matrices CR^.13 
and CC^.3 are unity matrix with a mult if i cation of a constant. 

Simplifying the above equation and putting the value of 
Cpi L2 +<C f /R * ' i L2' 5 ■ from equation (III. 3). 

PiL 1 + 

p i L4 + 

[ CB L14 3 Xl2 3CK 7 3R+CB L14 3TcC P CB L11 3 ]i L 1 + 

[ 

Above equation can be simplified by using equation (III.1) 

I 

pi. = LK..H, .+ CK 3i + CK j ,_ 3V I5j .+ CK.-Ue (IV. 4) i 

—LI 11 — L4 12 -LI 13 — B1 14 - 

Where, 

CKp= CB L14 3 T CB L12 3CK 6 3R+CB L14 3 T CR f 3CB L14 3 

“lO 3 " CB L14 3T,:B L12 3 ' :K 7 3R+[:B L 1 4 3T,:C P l:B L11 3 | 

CK 11 3 ^ ,:K 8 :I " 1 [' :B L14 3T|:B L12 3 ' :K 8 3R+ ' :B L14 3Tl:C P ,:B L14 3 ] j 

tK 12 3i-CK 8 D _1 fi:B L14 3 T i:B L12 DCK 5 DR+t:B L14 3 T CR f ]CB L14 3]cK 1 3+CK 10 3 i 



CK 13 3=-CK a 3- 


[ Ct >L 14 3T 


CBL 12 aCK 5 3R + C B L14 3 ' CR f 3CB L1 ,3 | CK^ 


'I 3 ] 


nK 14 D=-CKp 1 JcB L14 n T CB L12 3CK 5 3R+CB L14 D T i:R f 3CB L1 ^ | CKgll 
From equation (3.24), 

P -B1 = CR f" ,p iB1 + CC f 3 -B1 






=CR f 3CB L11 3pi L1 +CB L12 3RCpi L2 +(C f /R)i L2 3+CR f lCB L14 3pJ. L4 + 


CCpCB. .Hi, .+CCIDCB. . A li , . 
f L11 —LI f L14 — L4 


Putting the value of Cpi ^ ? +(C^/R)i^^3 from equation (111.3) 
P^B1 = -[( C;R f :,,:B L1i : " > iL1 +tB L1Z :,R [ i:K 5 :1 PiL1 +,:,< 6 : " > iL4 +,:K 7 :, i 


LI 


CK 8 :, iL4 + 


•]«R # 3C B Li 4 3pi L4 *CCJ3i: B L1 i31 L1+ CCJaCB Ll4 31 L4 ]] 


[‘ 


CR f lCB L11 3+CB L12 DCK 5 3R 


l] pi L1 -[cR f] CB L „3 + CB L1E DCK 6 3 R ] pi^- 


[ ™L 1E 3CK 7 3 R «cpCB Ll1 3]l L1 -[c B Li 2 3n< 8 3R*t c J3CB Ll 43]lL 4 


S£> 


P^B1 =i— f BK 19^P'ii_1 +BR E0^ P -L4 + ^ K E1^— L1 +BR 22^— L 4 ] <III ' 5> 

Where, 

CK 19 3=CR f 3CB L11 l+CB L12 DCK 5 3R 

CK 20 3=CR f 3CB L14 3+CB L12 lHK 6 DR 



Using equation (III.1) and (III.4),the equation can be simplified 
as follows: 


+[nK 19 ccK 2 ]+i;K E0 ]i:K 13 D]y B1 +[i:K 19 3CK 3 3+t:K zo 3CK 14 aJeJ 

or r 

pV _ * -CK.-Ui. -CK.-DV- -CK.-He (III. 6) 

— B1 15 —LI 16 — L4 17 — B1 18 - 

Where , 

CK 15 3=CK 19 nCK 1 l+CK 2() lC:K 12 D+CK 21 3 
CK 16 1=CK 20 3CK 11 3+CK 22 3 
CK iy 3=CK 19 1CK 2 3+CK 20 3CK 13 3 
CK 1S 3=CK 19 3CK 3 3+CK 20 3CK 14 3 



APPENDIX IV 


= 94.11 Volts line to ground 
= 0.158 ohms per phase 
= 0.073 mH per phase 


The datas of source, converter and dc load for digital 
simulation of the GTO converter are given below s 
(a) AC source data 

(i) Ac system frequency = 50 Hz 

(ii) Source voltage 

(iii) Source resistance 

(iv) Source inductance 
Cb) Converter data 

Ci) Overvoltage suppresion circuit elements 
Resistance = 6 ohms 

Inductance = 8 jjF 

Cii) Number of pulses per cycle CP) = IE 
(iii) Modulation index for EPWM = 0.8 

(c) Load data 

(i) Smoothing resistance (Rd ) =3.1 ohms 

(ii) smoothing inductance (Ld) = 24.0 mH 

(iii) Back ep»f (Vc) = 160.0 volts 

These datas are for rectifier operation . For inverter 

j 

operation all other datas are same except back emf which is taken; 

I 

as -180.0 volts. 1 


1 



APPENDIX V 


MODE SEQUENCE OF RECTIFIER OPERATION FOR m=O.S AND P= 12 






INTERVAL 


PERIOD IN DEG. 


G_ G„ G . G„ G , 

2 3 4 5 6 



“10" /V io 
^ 10 " *11 
“ll“ ^11 
^ll” 01 12 
" 12 " *12 
^12 _ 01 13 

-i- /» 1 


* 24 degree indicates the power mode of converter 
## tP degree indicates the free wheeling mode. 

X shows the GTOs conduction for a given period. 
They are related by following equation ' 


POWER MODE + FREE WHEELING MODE 







APPENDIX VL 


SYSTEM DATA (adapted from reference [ 14 ] ) 


DC base voltage = 1 00 kV 

DC base current = 1 kA 

DC base impedance = 100 Ohms 

Rated Power of DC Link = 240 MW for 12 pulse system 

= 120 MW for 6 pulse system 

AC system frequency = 50 Hz 


Bridge Transformers 

Resistance R . = R ^ 0 = 0.5 Ohms 

cl c z 

Commutating Reactance X c1 = X c2 


6.283 Ohms 


Smoothing Reactors 

Resistance R-. = R = 0.1 Ohms 

01 02 

Reactance X d1 = = 314.159 Ohms 


Transmission Lines 
Total resistance = 8.64 Ohms 

Total Inductance = 0.5014S H 

Total capacitance = 54.1645 pF 

Any number of PI sections upto a maximum of 10 may be specified 


Digital Controller (adapted from reference ) 

a) Constant Current Controller 

= O.^ radians/per unit current 


b) 


"1 

o 

, 2 = O.j^ radians /per unit current 


INVERTER CONSTANT C OWfcHT CONTR 
o '> P-12. 

Kl ~ • COHSTAVYT 0,0 

Tj =• 0.025 > 



•f- iu u ;.'Ui 

6 2 1 .33 1 5 325.^iit>S90 L 

$A,£l*F Ms book is to be returned on the 
date last stamped. 
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